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We present a first-principles computational scheme for investigating the ballistic transport proper- 
ties of one-dimensional nanostructures with non-collinear magnetic order. The electronic structure 
is obtained within density functional theory as implemented in the full-potential linearized aug- 
mented plane-wave (FLAPW) method and mapped to a tight-binding like transport Hamiltonian 
via non-collinear Wannier functions. The conductance is then computed based on the Landauer 
formula using the Green's function method. As a first application we study the conductance be- 
tween two ferromagnetic Co monowires terminated by single Mn apex atoms as a function of Mn-Mn 
separation. We vary the Mn-Mn separation from the contact (about 2.5 to 5 A) to the far tunneling 
regime (5 to 10 A). The magnetization direction of the Co electrodes is chosen either in parallel 
or antiparallel alignment and we allow for different spin configurations of the two Mn spins. In 
the tunneling and into the contact regime the conductance is dominated by 8-0(22 -states. In the 
close contact regime (below 3.5 A) there is an additional contribution for a parallel magnetization 
alignment from the dxz- and dyz-states which give rise to an increase of the magnetoresistance as 
it is absent for antiparallel magnetization. If we allow the Mn spins to relax a non-collinear spin 
state is formed close to contact due to the competition of ferromagnetic coupling between Mn and 
Co and antiferromagnetic coupling between the Mn spins. We demonstrate that the transition from 
a coUinear to such a non-collinear spin structure as the two Mn atoms approach leaves a character- 
istic fingerprint in the distance-dependent conductance and magnetoresistance of the junction. We 
explain the effect of the non-collinear spin state on the conductance based on the spin-dependent 
hybridization between the d^^^^y^-states of the Mn spins and their coupling to the Co electrodes. 



I. INTRODUCTION 

Break junction experiments have allowed to perform 
transport studies on nanoscale metallic contacts in which 
the mean free path of the electrons is much larger than 
the junction length. The observation of quantized con- 
ductance in such systems is a hallmark of ballistic trans- 
port and opened new vistas to study the scaling of elec- 
tronic devices down to the atomic length scale.i^ A draw- 
back of such experiments is the limited control of the 
microscopic arrangement in the junction which hinders a 
straight forward interpretation of the data and makes a 
comparison with theoretical calculations difficult.^ In this 
respect, a great advantage is given by the use of scanning 
tunneling microscopy (STM) experiments, in which a tip 
can approach and contact single atoms or molecules on 
a surface.l^^ In such experiments, it has been possible to 
measure the conductance as a function of tip-sample dis- 
tance from the tunneling to the contact regime. Due to 
the promise of spintronic devices for future applications 
with low power consumption and high speed, a recent 
focus of such contact measurements has been magnetic 
systems, e.g. spin-valve behavior has been obs erve d in 
single magnetic molecules or atoms on surfaced^ and 
the occurrence of the Kondo effect has been found in fer- 
romagnetic atomic contacts. 

It has been emphasized that the low coordination 
of the contact atoms in nanoscale junctions leads to 
an enhanced tendency towards magnetism, e.g. mag- 



netic moments are formed in systems of otherwise non- 
magnetic materials .'^^^^ Naturally, transport phenom- 
ena in such magnetic low-dimensional systems have 
raised a lot of attention and triggered many theoret- 
ical studies, which mainly focused on systems with 
coUinear magnetic order, considering also the effect 
of magnetoresistance.'iSJEi] \i .^^^s also recently realized 
that, if the magnetization direction of the two elec- 
trodes is opposite, a domain wall can form in the 
contact between them and the non-collinear order in 
the domain strongly affects the conductance and the 
magnetoresistance.'^2H2l] Finally, the effect of spin- orbit 
coupling on the conductance needs to be considerecP^I^Il 
which leads to novel transport pheno mena such as the 
ballistic anisotropic magnetoresistanc^^^ or the tun- 
neling anisotropic magnetoresistance.'^^^'^Sl 

Recently, the transition regime from tunneling to con- 
tact in a spin-polarized STM geometry has been stud- 
ied based on density functional theory in order to ex- 
plain e.g. the conductance of a single magnetic atom,'^ 
and to analyze the contribution from different conduc- 
tion channels.!^ As a magnetic STM tip approaches a 
single magnetic atom on a surface an exchange interac- 
tion with the tip apex atom occurs. In principle, it is 
possible to switch the magnetic moment of the adatom 
in such a way.l^^ If the magnetic moment of the adatom is 
exchange coupled to the substrate (as in Ref. |8j) there is 
a competition of exchange interactions which can result 
in a canting of the spins close to contact. Non-collinear 
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spin alignment in such an atomic contact can also occur 
if the adatom spin is canted due to exchange coupling on 
a substrate with a spin spiral structure as in Refs.[5Uland 
[3TJ The effect of such a non-collinearity in the spin direc- 
tion of the tip apex and the adatom on the conductance 
is the focus of the present work. 

We introduce an approach to calculate the conduc- 
tance in magnetic nanoj unctions with non-collinear spin 
structure from first-principles, employing the methodol- 
ogy of non-collinear Wannier functions (WFs) , which we 
describe in detail. In order to start from an accurate 
description of the electronic and magnetic structure of 
the system we use the full-potential linearized augmented 
plane wave (FLAPW) method based on density func- 
tional theory. We map the electronic structure of a sys- 
tem in a non-collinear magnetic state from the FLAPW 
description to a tight-binding like Hamiltonian via WFs. 
Finally, we calculate the conductance within the Lan- 
dauer approach with the technique of Green's functions. 

As a model system, we consider two Co monowires to 
each of which a single apex Mn atom is attached. We vary 
the distance between the two Mn atoms in order to cal- 
culate the conductance from the tunneling to the contact 
regime. The magnetization direction of the two Co elec- 
trodes is chosen either parallel (P) or antiparallel (AP) 
which allows us to obtain the distance-dependent magne- 
toresistance. In the tunneling regime the conductance is 
dominated by states of s-d22 -orbital character and only 
in the contact regime there is an additional contribution 
due to (ia^z^yz-states. As the latter conduction channel is 
suppressed in the AP alignment the magnetoresistance 
displays a large rise close to contact. 

When the two Mn atoms approach in the P electrode 
alignment a competition of the exchange interactions be- 
tween the two Mn spins and the Mn spins with the Co 
electrodes occurs. While the Mn spins couple ferromag- 
netically to the Co electrodes, they couple antiferromag- 
netically with each other. As a result a non-collinear 
arrangement becomes the magnetic ground state and the 
angle between the two Mn spins changes gradually from 
zero to about 105° at the closest separation we consid- 
ered. The conductance displays a characteristic dip as 
the non-collinear state forms which is also apparent in the 
distance-dependent magnetoresistance. We explain this 
reduction of the conductance due to non-collinear spin 
states from the spin-dependent hybridization of d^^.yz- 
states between the two Mn atoms which depends on the 
angle between their spin moments and partly suppresses 
the conduction in this channel. 

The paper is organized as follows. In Sec. |lT] we in- 
troduce our method to calculate the conductance of a 
one-dimensional nanoscale junction with a non-collinear 
spin structure. We discuss the extension of Wannicr func- 
tions to systems with non-collinear order (Sec. II A), the 
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implementation within the FLAPW method (Sec. 
and the incorporation into our transport code (Sec 
In Sec. |III| we introduce our model system consisting of 
two Co monowires to each of which a single Mn atom is 



attached. First, we analyze the magnetic and transport 
properties of collinear spin states from tunneling to con- 
tact (Sec. Ill B ) before we address the occurrence ofnon- 
collinear spin states in the contact regime (Sec. IIICl. 
We analyze the ballistic conductance of such spin states 
(Sec. Ill D I and show that a characteristic fingerprint is 
observed in the distance-dependent conductance and the 



magnetoresistance (Sec. HIE I. We end with a summary 
in Sec. HVl 



II. METHOD 

The density functional theory (DFT^ states that 
the energy functional of a general magnetic system 
-E[n(r), m(r)] is uniquely determined by the charge den- 
sity n(r) and the magnetization density m(r). The most 
common approximation made to a general magnetic sys- 
tem is to assume a collinear magnetization density, i.e., 
m(r) — r7i(r)e, where e is an arbitrary direction. Within 
this collinear approximation the energy is a unique func- 
tional of the charge density n(r) and the scalar magne- 
tization density m(r). Due to decoupled spin and real 
space, the spin-channels can be treated independently. 
However, it is known that relaxing the collinear approxi- 
mation and allowing for non-collinearity of the magneti- 
zation density in real space in the DFT setup leads to an 
ability of reliably treating whole classes of new phenom- 
ena, which rely on the properties of complex magnetic 
states.l^Sl 



A. Non-collinear Wannier Functions 

Within the DFT formulation for non-collinear mag- 
netic systems one solves the Kohn-Sham equationJ22 
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km/ ' 



(1) 



where I2 is the 2x2 unity matrix, V is the potential 
matrix which also mixes the spin channels, m is the band 
index and = (l^kmt) ' IV'kmi))^ is the spinor Bloch 

function with spin- up and spin-down components iV'kmt) 
and \i>kmi), respectively. 

For the M converged spinor Kohn-Sham orbitals 
\'*Pkm) ^ uniform mesh of Af k-points, the orthonor- 
mal set of Wannier functions can be obtained via the 
transformatior|23 
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m — 1 



(2) 



where the number of WFs N is smaller than or equal 

(k) 

to M and the matrices Umn represent the gauge free- 
dom of the WFs. In case when N — M and the group 
of bands we are extracting the WFs from is isolated 

(k) 

from other bands, the J7mn-matrices are unitary at each 
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fc-point. Imposing the constraint of maximal localiza- 
tion of WFs in real space determines the set of Umn- 
matrices up to a common global phase, and the corre- 
sponding WFs are called maximally-localized Wannier 
functions (MLWFs).E3Like the Kohn-Sham orbitals, the 
MLWFs are spinors and can be written as |VFkm) = 
(|VFkmt) I |W^km^))'^ in terms of their spin-up and spin- 
down components |Wkmt) ^nd |W^km^), respectively. 

For the construction of MLWFs within DFT elec- 
tronic structure codes, the matrix elements Mmn''' — 

(V'kml e"*''''" IV'k+bn) and aI^I = (V'kmlffn) need to be 
computed, where \gn) is a localized orbital, which defines 
the starting point of the iterative procedure of determin- 
ing the MLWFs. Since spin-up and spin-down are cou- 
pled in noncollinear calculations, these matrix elements 
involve a summation over the spin a. 

The matrix elements H„„/ (Ri — R2) for the WFs tight- 
binding Hamiltonian 

HwFs = 5m H„„'(Ri - R2) |Wr,„) (Wr,„.| (3) 



are given by^SI 
H„„' (Ri — R2) 
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.(4) 

Even though the Wannier and Bloch functions arc spinor 
valued, the transformation of the Hamiltonian from 
Bloch into Wannier representation is fully determined by 
the matrices Umn and the eigenvalues CmO^) as in the 
coUinear case. 



B. Non-collinear Wannier Functions within the 
FLAPW method 

The treatment of noncollinear magnetism within the 
FLAPW metho d as implemented in the Jiihch DFT 
code FLEUfP^I^ neglects the effect of intra-atomic non- 
coUinearity. Space is partitioned into the muffin-tin (MT) 
and interstitial regions (IR) . The spin density m(r) in the 
IR is treated without shape approximation as a continu- 
ous vector field. In the MT-sphere MT" of atom a only 
the projection of the spin density onto the direction e^j 
of the average spin moment is used for the generation of 
the exchange-correlation potential. The explicit one- and 
two-dimensional implementations also contain a third re- 
gion, the vacuum region (VR), which can be treated anal- 
ogously to the lR| 38 | 39 lThus, the self-consistent spin den- 
sity is approximated as 



m(r) 



m(r) 

m"(r)e' 



M 



IR (VR) 
MT" 



(5) 



A part of the intra-atomic non-coUinearity can still be 
described within this hybrid approach by decreasing the 
MT radii. 

The radial solutions (r) of angular momentum I of 
the scalar-relativistic Schrodinger equation in MT" and 
their energy derivatives uf^a (r) are calculated for the two 
spins cr" and used for the expansion of basis functions and 
Bloch functions. The spin quantum number cr" refers to 
the local spin quantization axis e^. The expansion coef- 
ficients of the eigenspinors of the local spin quantization 
axis in terms of the eigenspinors of the global spin quan- 
tization axis, which is the z axis, are given by 



- (expH|)cos(f),exp (^|)sin(f))^ 



X"^ = (-exp(-if)sin(|),exp(if)cos(|))'^. 



(6) 



where (j) and 9 are azimuthal and polar angles of the spin 
direction of MT" with respect to the global frame g. 

Within the MTs the wave function Vkm(i") is thus 
given by: 
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(7) 

where L denotes the angular momentum quantum num- 
bers and m is the band index. 

Using functions \gn) which are restricted each to a sin- 
gle MT-sphere has been found to result in a very good 
starting point for the iterative optimization of collinear 
WFs. "^3 Due to the approximate intra-atomic coUinear- 
ity, it is reasonable to choose in the noncollinear case the 
localized orbitals |g„) to be eigenstates of the projection 
of the spin-operator onto the local spin-quantization axis 

5n(r)|MT°(") ^'^CnLUni{r)YLir))x'^l"]^ . (8) 

L 

Here, c„l are expansion coefficients, a(n) is the index of 
the atom for which gnir) is non-zero, cr(n) is the spin 
associated with this trial orbital (7„(r) and Uniir) is the 
radial part of the trial orbital. Thus we obtain: 



A. 



k) 



L 

X { i^L^in) (k)] * / "riS ir)u,u {ry dr (9) 



The MT contribution to the Mmn"^ matrix may be 
written as 

^i^n'^^lMT^ V / d3r[V'km.= (r)]*e-^b-"-Vk+b„."(r), 



(10) 
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with V'km<T° (r) given by 



L 



(11) 



X [A^^,. (k)<.„ (r) + B!^^,^ (k)^r.= (0)] . 



The computation of the Mmn''^ matrix for noncolhnear 
systems reduces therefore to integrals for which expUcit 
expressions have been given for the FLAPW methodP^ 



C. Ballistic transport in systems with non-collinear 
magnetism 

The extension of the coUinear scheme for balhstic 
transport, which we described in detail in our previ- 
ous publication, Ref. I16[ to a non-collinear setup is now 
rather straightforward. Given the minimal WFs Hamil- 
tonian Hwfs (Eq- ([sj) for a non-collinear system, we are 
able to construct the tight-binding Hamiltonian of the 
nanojunction in accordance to our transport method,!^ 
which employs the partitioning of space into the scatter- 
ing region (5), as well as left (L) and right {R) leads: 



H 




(12) 



Compared to a coUinear calculation we have to deal with 
twice as many Wannier functions due to the insepara- 
ble spin-channels. Two calculatio ns using the locking 
techniquelSl are required for Eq. (12 1. The Hamilto- 



nian matrix Hg and the matrices H^s, describing the 
coupling to the leads, are obtained from a super-cell 
FLAPW calculation. The Green's functions for the leads 



GL/R{E) = [{E + te)I 



L/R 



Hl/h] ^ can be brought to 



finite sized surface Green's functions Sl/r by construct- 
ing /ji based on principal layers h^, /r and the coupling 
matrices hi^i^ / fifi^ Those matrices are obtained from a 
separate calculation of a perfect periodic lead. Following 
our Landauer-Biittiker method the non-collinear ballistic 
transport can be calculated with the Green's function of 
the scattering region: 

G5(^) - [Els - Hs - UlggLtlLS - lil,,gRllsR]-\ 

(13) 

The interaction between scattering region and the leads, 
and the resulting level broadenings are described by 
broadening matrices F 

Tl/r{E) - i[^l/r{E) - S[/^(i?)], (14) 

where Si/^(i?) are the self-energies of the leads: 

IIl/rXE) = fll^^^^gL/RiE)ULS/SR- (15) 

Finally the ballistic transport process is described by the 
transmission function T{E) 

T{E) = T^[Gs{E)Tl{E)G\{E)Tr{E)]. (16) 



resulting in the conductance through the junction 

.2 1 



G{E) 



-T{E) = -G^T{E) 



(17) 



with the conductance quantum Go = 2e^ /h. In the non- 
collinear case the trace operation of Eq. ([16]) has to be 
additionally performed over the spin a. The spin-channel 
information is therefore lost for general non-collinear sys- 
tems. 

We tested non-collinear Wannier functions on free- 
standing non-collinear magnetic Mn chains and found 
them to reproduce the FLAPW electronic structure with 
any given accuracy. For the performed transport cal- 
culations a 2"^^ nearest-neighbor (NN) tight-binding like 
Hamiltonian is sufficient due to the excellent correspon- 
dence of FLAPW and WFs electronic structure in the 
vicinity of the Fermi level in absence of s-d^^ band edges 
in that particular region, which usually would require to 
consider more neighborgi^. For the system, considered in 
the following, the orbitals participating in transport can 
be arranged according to the symmetry into the Ai (s 
and 1^22 orbitals), the A3 {d^z/yz) ^-nd the A4 {d^y/x'^-y'^) 
groups. 

III. THE Cooo-Mn-Mn-Cooo JUNCTION 

In the following section, we investigate the ballis- 
tic transport properties of coUinear and non-collinear 
magnetic configurations of a Cooo-Mn-Mn-Gooo junction, 
consisting of semi-infinite ferromagnetic Go monowires 
with magnetic Mn "tip" atoms, see sketch of the struc- 
ture in Fig. [1] We will discuss the effect of non-collinear 
magnetism on ballistic transport through such a junction 
specifically keeping in mind tunneling-to-contact STlVpWSl 
and mechanically-controllable break-j unction experi- 
ments. In particular, we investigate the changes in the 
transport properties upon changing the distance between 
the two Mn atoms, while keeping all other interatomic 
distances fixed at their equilibrium "semi-infinite" val- 
ues. We will show, that upon bringing the leads together 
the non-collinearity in this system emerges as a result of 
competing Mn-Mn and Mn-Go exchange interactions. 

We then explore the influence of non-collinear mag- 
netism on ballistic transport for various coUinear and 
non-collinear configurations. The nomenclature for the 
magnetic states in the Gooo-Mn-Mn-Gooo junction in- 
cludes the alignment of the magnetization directions of 
the leads, parallel (P) or anti-parallel (AP) to each other, 
and the directions of the two Mn spins. Without loss of 
generality, these directions are denoted with respect to 
the left lead which has magnetization "up". The Mn 
spins can point "up" (f), "down" (|) or in a direction 
which makes an angle a with the direction "up", see 
Fig. [1] In the latter case we consider the symmetric con- 
figuration, denoted as Pa, in which the spins of the Mn 
atoms make an angle of 2a between each other. For all 
considered non-collinear Pa states we fixed the direction 
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FIG. 1: (color online) Calculated magnetic configurations in 
the Cooo-Mn-Mn-Cooo junction. The magnetization of the left 
lead always points up (t) . The magnetization of the right lead 
can be either parallel (P) or anti-parallel (AP) to it. The di- 
rection of the Mn spins is marked with respect to the left lead: 
pointing up (f), down (J,), or along a direction at an angle a 
with the magnetization of the left lead. Special configurations 
of interest which are discussed in the text are shown: (a) Pfti 
(b) symmetric Pa state with a = 105°, P105°, (c) PU (d) 
Pit and (e) APfi. In the text, we refer to the separation 
between the leads in terms of the distance between the Mn 
atoms, dMn-Mn. Co atoms are displayed as orange spheres, 
Mn atoms are displayed as red spheres. 

of all Co atoms either up or down, depending on the 
magnetization direction of the corresponding lead. The 
energy differences between different magnetic states are 
given per Mn atom. 



A. Computational details 

For all collinear and non-collinear electronic struc- 
ture calculations we used density functional theory 
within generalized gradient approximation (GGA) to the 
exchange-correlation potentialj^Sl as implemented in the 
FLAPW Jiilich code FLEUR.'23 The wires were calculated 
in three-dimensional super-cells, with an interchain sepa- 
ration in the x-y— plane of 13 bohr. The super-cell setup 
along the chain's axes (z-direction) is described in detail 
below. The Brillouin zone (BZ) was sampled by 12 or 
24 fc-points along the z-axis, depending on the size of 
the super-cell. All calculations were performed with an 
LAPW basis cut-off parameter fcmax of 3.7 bohr~^, re- 
sulting in approximately 625 LAPW basis functions per 
atom. 

The parallel magnetic configuration (P) of Cooo-Mn- 



Mn-Cooo junctions was investigated in an 8-atom super- 
cell along the chain direction, consisting of six Co atoms 
with an equilibrium interatomic distance of the Co in- 
finite monowire of dco = 4.15 bohr, and two attached 
Mn atoms, see Fig. [l] For all considered magnetic con- 
figurations, with parallel or antiparallel alignment of the 
magnetization of the leads, as well as non-collinear mag- 
netic states, irrespective of the separation between the 
leads, we fixed the Co-Mn distance dco-Mn to 4.48 bohr, 
which corresponds to the equilibrium distance between 
the ferromagnetic Co and Mn atoms at a very large sep- 
aration between the leads. For the Pft and P|J,-states 
of the junction we considered the inter-Mn separation 
of dMn-Mn = 5.0, 5.5, 7.0, 8.5, 10.0, 12.5, 15.0, 17.5 
and 20.0 bohr. The anti-parallel magnetic configuration 
(AP) of Cooo-Mn-Mn-Cooo junctions was calculated in 
a 16-atom super-cell, consisting of six ^-Co atoms and 
two Mn atoms on one side, and six J,- Co atoms and two 
Mn atoms at another end of the junction. In this case 
^Mn-Mn was Set to 4.5, 5.0, 5.5, 7.0 and 8.5 bohr. 

For the conductance calculations we applied the lock- 
ing technique to a perfect monowire to describe the semi- 
infinite leads, as described in detail in Ref.[T6l In all cases 
the Wannier functions were generated on a 1 x 1 x 24 fc- 
point grid in the BZ. For the collinear cases the WFs were 
generated from one 4s- and five 3(i-orbitals per atom for 
each spin separately, which were constructed from the ra- 
dial solutions for the FLAPW potential. In non-collinear 
calculations the spin channels are mixed, and two 4s- 
and ten 3(i-orbitals per atom were used to construct the 
WFs per atom. The energy bands were disentangled us- 
ing the procedure described in Ref. 021 For the collinear 
calculations the lowest 58 eigenvalues per fc-point were 
used to obtain 48 WFs for the 8 atom super-cell and the 
lowest 104 eigenvalues per /c-point for 96 WFs for the 16 
atom super-cell calculations. With non-coUinearity of the 
magnetization included the lowest 103 eigenvalues per k- 
point were used to obtain 96 WFs for the 8 atom unit 
cell. For testing purposes, for several non-collinear con- 
figurations we compared the electronic structure of the 
system calculated with FLEUR and with corresponding 
WFs, finding that a very good description of the elec- 
tronic structure can be achieved with WFs within the 
3rd nearest-neighbor approximation, while for WFs cal- 
culations of the transmission in the vicinity of the Fermi 
level already the 2nd nearest-neighbor approximation to 
the WFs Hamiltonian provides very reliable results. 



B. Collinear magnetic states of the junction from 
tunneling to contact 

We start the investigation of the Cooo-Mn-Mn- 
Cooo junction with both leads positioned far away from 
each other. To mimic a tip-sample approach, we decrease 
the Mn-Mn distance rfMn-Mn, and calculate the energies 
of the collinear states Pfti Pii and P\.'t, showing the re- 
sults in Fig. [2];a). The energy difference between Ptt and 
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P^l states when the distance is varied in the tunnehng 
regime from divin-Mn = 20 bohr down to 10 bohr re- 
mains relatively constant and constitutes around 27 meV 
per atom, indicating weak interaction between both sides 
of the junction and a weak ferromagnetic coupling be- 
tween the Mn atom and its nearest Co neighbor (NN 
Co). After a small reduction of the energy difference be- 
tween the Pit and states around dMn-Mn = 8 bohr, 
the ferromagnetic (FM) Mn-Co coupling becomes more 
stable for decreasing dMn-Mn, expressed in an increas- 
ing energy difference. In the contact regime we find a 
slight decrease in the energy difference from 93.5 meV 
at dMn-Mn = 5.0 bohr down to 77 meV per atom at 
^Mn-Mn — 4.5 bohr. This decrease in energy can be cor- 
related with strong changes in the Mn and NN Co spin 
moments, ^Mn and S'mii, respectively, upon decreasing 
the distance, see Fig. ^h) and (c) (see also discussion in 
the next section). 

While in all cases the spin moments of the Co atoms, 
not neighboring the Mn atoms directly (sa 2.09 /is) are 
very similar to the spin moments of the Co atom in 
an infinite lead (w 2. 07 /is), the spin moments of Mn 
atoms and the NN Co atoms can be strongly affected by 
dMn-Mn at closc contact and the spin configuration of the 
junction. Namely, for the PJ,| state S'mii decreases from 
4.3 fiB to 3.9 fiB, while NN Sco increases from 1.4 
to 1.8 /is, as dMn-Mn is varied from 5.5 to 4.5 bohr. On 
the other hand, if the Mn spin moment exhibits a simi- 
lar variation as a function of distance for the Pft state, 
the spin moment of the NN Co atoms remains relatively 
constant (w 2 /is). This interplay between structure and 
magnetism already indicates that the intra-atomic as well 
as inter-atomic exchange, given by the Stoner parameter 
/ and the Heisenberg exchange constants J, respectively, 
may be of importance for further understanding of the 
magnetic properties of this system. 

The change from FM coupling at larger interatomic 
distances to an antifcrromagnetic (AFM) coupling at 
smaller dMn-Mn in an infinite Mn chain has been pre- 
viously predicted based on DFT calculations.^'^ In the 
vicinity of this crossover point the Mn spins favor non- 
coUinear magnetic order .1^31^ To demonstrate a strong 
tendency of Mn spin moments to AFM coupling at 
smaller values of dMn-Mn we plot the energy difference 
between the Pft and P|t states in Fig. ^a). Revers- 
ing one of the Mn spin moments in the Pft configu- 
ration is clearly energetically more favorable than the 
Ptt state when the distance between the Mn atoms is 
below « 5.2 bohr. In this case the gain in energy due to 
switch of the Mn spin moment can be explained only by 
the strong AFM coupling of the two Mn atoms for this 
regime of interatomic distance, since the coupling of the 
Mn atom with its NN Co atom is ferromagnetic. The 
spin moments in the P|t state at rfivin-Mn = 4.5 bohr 
constitute 4.1 /is for Mn and 1.5 /is for its NN Co on the 
AFM side, and 3.7 /is for Mn and 1.9 /x^ for the n.n Co on 
the FM side. For larger dMn-Mn values the Pft state is 
the lowest in energy as compared to all possible coUinear 
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FIG. 2: (color online) (a) Energy of the P4,4, state (red di- 
amonds, dashed line) and P^t state (green triangles, dotted 
line) with respect to the energy of the Ptt state, as a function 
of the Mn-Mn interatomic distance dMn-Mn. Spin moment of 
the Mn atoms, Smu, (b), and NN Co atoms, ^co, (c), for 
Ptt (black circles, solid line) and P^^ (red diamonds, dashed 
line) configurations are given as a function of dMn-Mn. (d) 
Total conductance at the Fermi level, Gtotai(-E'F) (upper line) 
and Ai conductance Gai(-Ef) (lower line) on a logarithmic 
scale for the Ptt-state. Gray shaded area is associated with 
the Ai conductance, while cyan shaded area with the A3 con- 
ductance. 



states of the junction in which the magnetization direc- 
tion of the left and right leads is the same, which is in- 
dicative of the FM Mn-Mn coupling for larger distances. 

In Fig. [2jd) we present the results of our calcula- 
tions for the evolution of the ballistic conductance of 
a Pft junction when going from the tunneling to the 
contact regime. The main contribution to the conduc- 
tance at large Mn-Mn distances is coming solely from 
the Ai channel, owing to the overlap between the s-d^2 
orbitals of the neighboring Mn atoms across the barrier. 
Within our approach, the expected exponential behav- 
ior of the conductance at very large distances is very 
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nicely reproduced. At a distance of divin-Mn ~ 10 bohr 
the conductance approaches the magnitude of the con- 
ductance quantum, reaching saturation upon further de- 
creasing the distance. For distances in the contact regime 
below 7 bohr more localized d-orbitals of A3 symmetry 
start contributing to the total conductance, as can be 
seen in Fig. ^d). The A3 contribution to the conduc- 
tance increases with decreasing distance. As we shall see 
in the following the details of hybridization between the 
A3 orbitals are very sensitive to the magnetic state of 
the junction. On the other hand, in all considered cases 
the d-states of A4 symmetry do not contribute to the 
conductance due to an energetic mismatch between the 
states of this symmetry of NN Co and Mn atoms, see 
discussion in section UlI Dl 



C. Non-collinear magnetic states of the junction in 
contact regime 



According to the findings presented above, we expect 
that Mn spin moments in the contact regime will experi- 
ence a frustration when the magnetizations of the leads 
are parallel to each other. In this case, when Mn atoms 
are close enough, FM coupling of Mn spins with NN Co 
atoms and AFM Mn-Mn coupling can possibly lead to a 
stable non-collinear magnetic state. In order to consider 
this situation, we introduce an angle a between the spin 
moments of the Mn and the NN Co atoms, rotating the 
first Mn spin moment by +a and the second one by —a, 
while keeping the moments of the Co atoms fixed, see 
Fig. [ijb). This is what we call a symmetric Pa-state. 
We choose a distance of 4.5 bohr between the Mn atoms 
as a representative of the contact regime at which the 
Mn-Mn coupling is strongly antiferromagnetic. 

The results of our calculations for the total energy 
of the Pa state, Epa, in relation to the energy of the 
Pft state are shown in Fig. [sjja) as a function of the an- 
gle a. From this plot we observe that the minimum of the 
total energy is acquired for the non-collinear P105° state, 
which is 137 meV lower in energy than the correspond- 
ing coUinear Pft state. The failure of a straightforward 
description of the energy landscape Epa in terms of a 
simple Heisenberg model which assumes just the nearest- 
neighbor Co-Mn and Mn-Mn exchange coupling, given by 
antiferromagnetic Jun < and ferromagnetic Jco > 0, 
respectively, can be understood from noticing that the ex- 
pression for the energy within this approximation, given 
by 



"(aT T 
400f«°-tD^ 



Epaia) 



1 



( Jmh cos(2a) + 2 Jco cos(a)), (18) 



acquires a minimum for angles a below 90°, in contra- 
diction to our calculations. 

The solution to this deficiency of the Heisenberg model 
can be given by lifting the assumption that the exchange 
interaction between the Mn and Co spins, given by Jcoj 
is ferromagnetic. As we can see from Fig. ^h) and (c). 




FIG. 3: (color online) (a) Total energy of the Pa state with re- 
spect to the Pft state as a function of angle a (solid line, filled 
circles). With open black squares we show the total energy of 
the Pa state from which the Stoner energy due to creation of 
the NN Co spin moments is subtracted. The corresponding 
fit of the corrected energy to the Heisenberg model is shown 
with a dashed blue line. With the red dotted line we show 
the total energy of the PtJ, state with only the ^-Mn spin ro- 
tated by an angle a, calcualted within the Heisenberg model 
and incorporating also the Stoner exchange, (b) Mn spin mo- 
ment as a function of angle a. (c) NN Co spin moment as a 
function of angle a. (d) Ballistic conductance at the Fermi 
level for the Pa state as a function of angle a: total conduc- 
tance GtotaiiEp) (upper line), the A3 conductance contribu- 
tion (cyan shaded area) and and Ai conductance Gai(£f) 
(lower line, gray shaded area). 



while the Mn spin moment remains relatively constant 
upon changing a, Sco for values of a below 60° is by as 
much as 0.2 /x^ larger than for a > 90°. Owing to the 
intra-atomic Stoner exchange, the non-collinear states 
with small a therefore acquire a negative contribution 
to the total energy in addition to that proportional to 
Jco, as compared to larger angles. If we account for 
the energy gain due to creation of the NN Co spin mo- 
ments by a Stoner parameter of Co, / « 990 meV,l^ and 
subtract the energy gain i^stoncr = '^^^Co from the cal- 
culated DFT dispersion, we arrive at the energy disper- 
sion (squares in Fig. [sjja)), which reflects only exchange 
interactions between the atoms. If we fit this curve ac- 
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cording to Eq.([l8]) (dashed line in Fig. [sja)), we obtain 
the " non-renormalized" Heisenberg exchange constants 
of Jqo = —170 meV and Jmh — —366 meV. It becomes 
clear now, that, although the "pure" exchange coupling 
between the Mn and Co spins is expectedly antiferromag- 
netic, the larger spin moment of Co in the parallel spin 
alignment with Mn, tips the balance in favor of ferromag- 
netic coupling between the spins, which can be observed 
for a large range of distances divin-Mn, c.f. Fig. [2]ja). 

In Fig. [2]ja) we also observe that, judging from the 
energies, in the close contact regime the collinear P^t" 
state is competing with the non-collinear Pa-state for 
the global ground state of the system. Indeed, our cal- 
culations show that at the dMn-Mn of 4.5 bohr the Pit- 
configuration is by a tiny value of 5 meV lower in energy 
than the P105° solution. We argue, however, that the 
P|t-state is not very likely to appear in experiments, 
given that the Co electrodes are identical. In this case, 
the adiabatic rise of the intrinsically asymmetric Pit- 
configuration via symmetric non-collinear states cannot 
happen, as the electrodes, initially being in the Ptt-state 
when very far from each other, are brought together (see 
also discussion at the beginning of section III E[ ). Nev- 
ertheless, it seems plausible, that such state, if observed 
in experiment, is created via a rapid flip of one of the 
Mn atoms in the contact regime, during, e.g., a reforma- 
tion of the lead geometry, or an inelastic current-induced 
spin-flip process. Our calculations, shown with a dotted 
line in Fig.|3](a), based on the Heisenberg model extended 
by the Stoner term of intra-atomic exchange of the Co 
moments, indicate, that once the system enters the Pit- 
state, it is effectively "trapped" there, since the i-Mn is 
energetically quite stable versus deviations in the angle 
its spin makes with the rest of the spins in the system. 
Thus, we do not consider any non-collinear states asso- 
ciated with the Pit-state in the following. 



D. Ballistic conductance of non-collinear magnetic 
states of the junction 

In this section we perform a detailed analysis of the 
balhstic conductance G{Ep) of the Pa state at fixed dis- 
tance between Mn atoms of 4.5 bohr. At this distance, 
we calculate G{Ep) as a function of angle a and present 
the results in Fig. |3jd). In this plot we observe that 
the conductance exhibits a very non-trivial dependence 
on a, originating mainly from the Aa-orbitals idxz,yz), 
while the Ai contribution (s — d^^) to the conductance, 
Gai, remains almost perfectly constant. Surprisingly, 
the Aa-conductance almost vanishes for a of about 70° , 
away from any high-symmetry spin configuration in the 
junction, suggesting, that the dependence of the details 
of hybridization and electronic structure on the angle be- 
tween the Mn spins can be rather delicate. In order to 
analyze this dependence in more detail, as a function of 
a, we plot the energy-dependent conductance, G{E), ver- 
sus the local densities of states (LDOS) of Mn and NN 



Co atoms resolved into spin-up and spin-down contribu- 
tions with respect to the global spin quantization z-axis. 
Fig. |4] Mainly we focus on the Aa-contribution to the 
conductance and the LDOS, and only in the upper (Ptt) 
and lower (Pii) panels of Fig. [i] we show also the total 
LDOS of the atoms. 

The conductance at a given energy E depends on the 
presence of available states in the LDOS of the atoms at 
E, and on the coupling between these states across the 
junction — both of which depend on the orientation of 
the spins with respect to each other. By looking at the 
LDOS of the atoms presented in Fig. |4] for a ~ 0° we 
can explain the absence of the A4 (dxy^x^-y^) contribu- 
tion to the conductance: the localized A4 states of the 
Co atoms, which can be seen as pronounced peaks in the 
LDOS marked with the dashed line in Fig. |4j are posi- 
tioned at about —2 eV for spin- up channel and directly 
at the Fermi energy for spin-down channel, while the cor- 
responding Mn A4 states are positioned below —2.5 and 
above -1-1 eV, prohibiting thus the hybridization between 
the Co and Mn orbitals of A4 symmetry across the junc- 
tion. Noticeably, the LDOS of both atoms for the up-spin 
in a wide region of energies around Ep is absent, lead- 
ing to a negligible t- conductance. Here, it is important 
to remark, that the LDOS of the NN Co atoms around 
the Fermi energy overall resembles quite well the LDOS 
of a Co atom in a Co monowire, see e.g. Fig. 13(a) in 
Ref. [ini or even of a Co atom deposited on a noble-metal 
surfaces, see e.g. Ref. US] This means that our results 
should be rather stable with respect to the geometry of 
the Co leads, manifesting that the main influence on the 
A3 conductance at Ep would come from the hybridiza- 
tion of the Mn and NN Co states. 

Turning now to the comparatively delocalized A3- 
states (solid line) on both Co and Mn atoms, we observe 
for Ptt (upper panels of Fig. |4| that they hybridize di- 
rectly at the Fermi energy, which leads to a significant A3 
contribution to the conductance. Specifically, while the 
A3 subband of Co spreads from —1.8 to -|-1 eV, the A3- 
down states of Mn atoms are distinctly split into wide 
bonding ("b") states at the Fermi energy and narrow 
anti-bonding ("a") states at -1-1.8 eV. Very importantly 
for the transport properties of the system, the hybridiza- 
tion of the A3 band of Co with the A3 states of Mn is 
non-trivial, (i) The A3 states of Co exhibit a dip at the 
position of the maximal density of bonding states of Mn 
due to the fact that these Mn states are localized mainly 
in between the Mn atoms prohibiting strong overlap with 
the Co states, (ii) The upper, antibonding, part of the 
Co A3 band hybridizes stronger with the bonding states 
of Mn, since the antibonding states of Co atoms have a 
larger overlap with the Mn orbitals, which results in a 
larger A3-conductance above Ep. (iii) Analogously, for 
energies below Ep the conductance is suppressed, since 
the bonding-like A3 Co states have smaller overlap with 
the Mn bonding states. 

Let us now follow the evolution of the electronic struc- 
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E-Ep (eV) E-Ep (eV) E-Ep (eV) 

FIG. 4: (color online) Transport properties and electronic structure of the Pa-state at dMn-Mn of 4.5 bohr as a function of 
angle a (indicated on the left), (a) Total conductance Gtotai (upper line) decomposed into the A3- (cyan shaded area) and 
Ai-contributions (Gai, lower line, gray shaded area). The spin-resolved local density of states (LDOS) of the NN Co atom 
and Mn atom are given in columns (b) and (c), respectively. For both (b) and (c) the LDOS spin-decomposition is performed 
with respect to the global frame, with spin-up and spin-down LDOS presented in the upper and lower parts of each plot. In 
both (b) and (c), the As-contribution is indicated with solid lines, while for a = 0° and a = 180° also the total LDOS is shown 
with dashed lines. The red triangles in (b) follow the development of the dip in the LDOS of the NN Co atom as the angle a 
is varied. In (c), the bonding and anti-bonding unoccupied Mn states are marked with "b" and "a", respectively. For details 
see text. 
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ture upon increasing the angle between the Mn spins. 
Two trends in the LDOS can be clearly observed in Fig.|4] 
Firstly, with increasing a the splitting between the bond- 
ing and antibonding Mn states decreases owing to the 
mixed spin character of the states. At the angle of 90°, 
when Mn spins are anti-parallel to each other, both types 
of states transform into degenerate As-orbitals of the 
"isolated" Mn atoms, since the hybridization between the 
Mn states of the same spin is almost absent due to large 
exchange splitting. On the other hand, the dip in the 
l-LDOS of the NN Co atoms follows the position of the 
bonding state of the Mn dimer, moving from the Fermi 
energy at a = 0° to -1-0.2 eV for a = 90° (indicated by 
filled triangles in Fig. |4| . Overall, such redistribution of 
the LDOS of the atoms combined with the effect of de- 
creasing LDOS of Mn atoms for spin-down channel at the 
Fermi energy when the angle a is varied, results first in a 
decrease of the conductance at Ep for a ~ 70°, followed 
by a consequent increase with increasing angle. 

When the angle a increases further beyond 90°, the 
bonding and anti-bonding Mn states eventually acquire 
their initial splitting bX a — 180° (P4,4,-state) , when the 
Mn spins are coUinear again. Simultaneously, with in- 
creasing angle, we observe that the Mn states around 
the Fermi energy become sharper, since the hybridization 
with the Co leads decreases as the Mn states become pre- 
dominantly spin-up in character. Interestingly, while for 
90° < a < 120° a large value of the Aa-conductance is 
due to a significant amount of delocalized Co and bond- 
ing Mn states at the Fermi energy in the spin-down chan- 
nel, for larger angles the value of Gi^.^(Ef) is due to a 
sharp resonant Co state in the spin-up channel at the 
Fermi energy, coupled to a bonding Mn state at Ep. 
When further increasing a above 170°, this resonance 
becomes more localized and decoupled from the states in 
the leads, while the Mn LDOS at the Ep in the minority 
spin-channel vanishes, causing a sharp drop in the A3- 
conductance. By looking at the total LDOS of the NN Co 
atom in the P^-state we observe that it remains basically 
unaffected, as compared to the Ptf-configuration, while 
the Mn states become pronouncedly decoupled from the 
states of the NN Co owing to the energetical mismatch 
for both spin channels. 



E. Fingerprints of non-collinear magnetic states of 
the junction in ballistic conductance experiments 

Finally, we investigate the evolution of the conductance 
and the magnetoresistance of different magnetic states of 
the junction within the contact regime mimicking a typ- 
ical STM or break junction experiment. Here, we are 
partly motivated by the fact that a non-trivial behavior 
of magnetoresistance when going from tunelling to con- 
tact has been recently observed in STM experiments, see 
e.g. Ref. m At a very large separation between the leads 
(or, the tip and the sample in the STM language), owing 
to the FM coupling of the Mn atom to the Co chain, one 



can imagine only two possible magnetic configurations — 
Pft and APfi. The conductance of these two magnetic 
states in the tunneling regime, arising mainly from the 
s-orbitals, is orders of magnitude smaller than in the con- 
tact regime, for which the dependence of G(Ep) on the 
distance can be non-trivial due to the large contribution 
of the c?-states. 

In the case of the AP^i configuration, the starting 
coUinear arrangement of the spins will survive over the 
whole range of the separation between the leads, since in 
the contact regime, when the Mn atoms are close to each 
other, both exchange preferences of the Mn spins, that 
is, FM coupling to the NN Co spins and AFM coupling 
among each other, are fulfilled. Small possible deviations 
from the coUinear arrangement of the Mn spins, which 
can affect the details of the distribution of the Aa-statcs 
and their coupling to the leads, would not manifest in 
a conductance measurement, owing to the antiparallel 
magnetizations of the leads, and corresponding complete 
dominance of the Ai channel for conductance at Ep in 
this case. Fig. [5] As we can see from this figure, G{Ep) 
lies in between 0.5 and 1.0 Gq, when the distance between 
the Mn atoms is varied from 8.5 to 4.5 bohr. This is very 
similar to the behavior of the conductance at the Fermi 
energy of pure AP Co leads without Mn atoms, see e.g. 
Fig. 11(b) of Ref.HSl 

Owing to the magnetic frustration of the Mn spins of 
the junction in the contact regime, for the Pft initial con- 
figuration, we consider the Pa and P|t states in addition 
to the Pft state when the Mn-Mn distance is relatively 
small. Here, as we have seen in the preceding section, the 
conductance at the Fermi energy can be very strongly in- 
fiuenced by the details of hybridization between the A3- 
orbitals. On the other hand, since very often transport 
measurements serve as the only experimental insight into 
the magnetic structure of a system, it is very important 
to coin each of the possible magnetic states with a unique 
fingerprint which can be related to the experimental data. 
Below, we suggest that indeed three distinct spin states 
in a Cooo-Mn-Mn-Cooo junction — which can occur in an 
experiment due to various reasons such as structural de- 
tails, temperature fluctuations, external magnetic field, 
etc. — lead to different transport signatures. 

As already shown in Fig. [2jd), the conductance of the 
coUinear Pff state rapidly rises towards a value of 1.8 Gq 
as the distance between the leads is decreased. Com- 
pared to other possible magnetic configurations of the 
junction, Gp^-f{Ep) is significantly larger in value, see 
Fig. [5] because of the alignment of the minority spin s- 
and d-states of the Co electrodes and the Mn atoms at the 
Fermi energy, which ideally favors perfect transmission. 
In contrast, the conductance of the coUinear P4,t-state is 
significantly suppressed, reaching only 1.0 Gq at the sep- 
aration of 4.5 bohr, due to the large exchange splitting of 
the A3 states of the Mn atoms with antiparallel spin mo- 
ments which hinders the A3 conductance. The conduc- 
tance of the non-collinear ground-state Pa-state lies in 
between the values for both limiting coUinear configura- 
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FIG. 5: (color online) Conductance at the Fermi energy of var- 
ious possible magnetic states of the Cocx)-Mn-Mn-Cooo junc- 
tion as a function of the distance between the Mn atoms. 
Following magnetic states are considered: APJ,')' (filled cir- 
cles, solid line), Pft (open circles, dashed line), P^t (dia- 
monds, dash-dash-dotted line), and Pa (squares, dot-dashed 
line). The Ai contribution for the P^t-state is shown with 
a thin solid line and grey shaded area, while the A3 part is 
shaded in cyan. For the Pa states the state which is lowest 
in energy among all possible angles a at a fixed distance is 
considered. In the inset the values of the magnetoresistance 
for different P-states are shown. 



tions. In the close contact regime, at dMn-Mn of 4.5 bohr, 
the conductance of the Pa-state of 1.4 Go is exactly in 
between the values of G'p^^{Ep) and G^\^^{Ef). Clearly, 
the difference of 0.4 Go, stemming from the variation in 
the As-conductance vi^ith the spin state, can be easily 
detected in experiment, allowing for a way to distinguish 
between different possible magnetic configurations. At 
the distance of 5.0 bohr the ground state among the Pa 
states is the P90° state, while at larger distances above 
5.5 bohr the system converges to a coUinear configura- 
tion. The angle a in the lowest in energy Pa state de- 
creases smoothly with increasing the separation, and we 
speculate, that owing to the non-monotonous behavior of 
the conductance as a function of a, seen in Fig. [sj^d), the 
conductance as a function of dMn-Mn, can exhibit several 
features similar to that at dMn-Mn of 5.0 bohr, although 
we did not perform the calculations to support this state- 
ment owing to the required computational effort. 

According to recent experiments)^ the conductance of 
the junction with the parallel (P) and antiparallel (AP) 
orientation of the lead's magnetization can be related to 
each other via measuring the magnetoresistance (MR). 
From the values presented in Fig. [5] we calculate the MR 



of the junction, defined as: 

MR = GAEf)-G^AEf) ^ ^QQ^^^ (^g^ 
Gap(£'f) 

and present the MR as a function of separation be- 
tween the electrodes in the inset of Fig. [5) where we 
choose Gap^^ Gap{Ef), and values of Gp^^{Ef), 
Gpi^{Ef) and GpaiEp) for Gp{Ep). The overall smaller 
AP4,t conductance as compared the P-configurations re- 
sults in positive magnetoresistance values. The MR 
curves as a function of the distance generally resemble 
those of the conductance, with the values of the MR of 
22, 62 and 105% at the distance of 4.5 bohr for the P|t, 
Pa and Pff-states, respectively. Much more pronounced 
in the MR is the feature characteristic to the V\.'\ and Pa- 
configurations — a dip around the distance of 5.0 bohr, 
also present in the conductance curves. As can be seen 
from Fig. |5] at this distance, the MR almost completely 
vanishes when the Mn spins exhibit a different from FM 
configuration. Overall, we conclude, that the pronounced 
difference in the shape and magnitude of the MR curves 
can be also used in experiments to shed light onto the 
complex magnetism in this type of systems. 

IV. SUMMARY 

In this work we presented the realization of a first 
principles scheme for calculating the ballistic transport 
properties of magnetically complex one-dimensional sys- 
tems employing the technique of non-coUinear Wannier 
functions. We use the FLAPW method in order to cal- 
culate the electronic structure of the system with high 
accuracy and use the Wannier functions to transfer it 
to our transport calculations performed within the Lan- 
dauer approach. As spin-orbit interaction can be nat- 
urally included into the consideration within this tech- 
nique, c.f. Ref. 1161 the method introduced here can be 
used to explore the rich field of transport phenomena in 
systems such as nano- or atomic-sized contacts, break 
junctions, or STM experiments for which both effects, 
spin-orbit coupling and frustrated exchange interactions, 
can be prominent. 

As a first application of our approach, we consider 
the ballistic transport properties of a single-atom junc- 
tion formed by two semi-infinite Co electrodes with a 
single apex Mn atom. We study the conductance as a 
function of the separation between the two Mn atoms 
from the tunneling to the contact regime taking into ac- 
count the complex magnetic interaction in the junction. 
As we demonstrate, even such a simple setup allows to 
draw some general conclusions concerning the interplay 
of structure and magnetism for the transport through 
such atomic-sized contacts which are in the focus of to- 
day's research. We analyze the ballistic conductance of 
the junction with lead magnetizations in parallel and an- 
tiparallel alignment. We consider separately the tunnel- 
ing (separation larger than about 5 A) and contact (below 
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5 A) regimes of the junction, and we demonstrate that in 
the tunnehng regime the conductance G is solely com- 
ing from the overlap between the Ai (s-(i22-orbitals) of 
the contacts. In this case the Mn spins prefer to order 
ferromagnetically with respect to the magnetization of 
the leads. On the other hand, upon reaching the close 
contact regime (below 3.5 A), when the magnitude of the 
conductance reaches 1 Gg, the hybridization between the 
A3 (ci^^z, (ij^^-orbitals) states of the junction starts to pro- 
vide a sizable contribution to G. 

In the close contact regime, when the hybridization be- 
tween the Mn atoms is significant, Mn spins experience 
a frustration due to the FM coupling with the leads and 
an AFM Mn-Mn coupling. The competition between the 
two gives rise to a stable non-coUinear solution which 
can be characterized by a tilting angle of the spins, a. 
General for this type of junction is the sensitivity of the 
d-orbital conductance on the angle a, which is due to 
a delicate interplay between the hybridization details of 
the Mn and Co states at the Fermi energy, as well as 
spin-asymmetry in their distribution. This gives rise to 
a non-trivial a-dependence of the conductance of the d- 
states. We show that the complicated Aa-channel con- 
ductance arising on the background of almost constant 
Ai contribution can be used in order to distinguish be- 
tween different magnetic states of the contact via either 
a direct conductance measurement, or via measuring the 
magnetoresistance, which, according to our calculations, 
can vary in the contact regime between 20 and 100%, 
depending on the spin arrangement. 

Finally, we would like to comment on our approxima- 



tion for the geometry of the junction we have assumed in 
this work. Albeit being very simple, it allows to capture 
the key features which govern the transport properties 
of the system, while keeping the computational burden 
reasonable. Namely, within this geometry: (i) the transi- 
tion from tunneling to contact can be naturally studied; 
(ii) the magnetic frustration of the spins in the junction, 
and (iii) the delicate details of the hybridization of the 
adatom with the lead reservoirs are taken into account; 
(iv) the sensitive dependence of the spin moments on the 
magnetic configuration in the nano-contact is included 
into our considerations. Of course, in order to achieve 
a quantitative agreement of the calculated values to the 
experimentally measured ones in this type of junction 
beyond the major trends, all details of the structure and 
structural reformation upon approaching should be ide- 
ally accounted for. Such a challenging study lies, how- 
ever, outside of the scope of the current work, and we 
leave it for future studies. 
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